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PACS. 1 - 64.70.Md Transitions in liquid crystals. 

PACS. 2 - 64.60. Qb Nucleation. 

PACS. 3 - 61.30.Vx Polymer liquid crystals. 



Abstract. - Single crystal nematic elastomers undergo a transition from a strongly ordered 
phase N to an "isotropic" phase /. We show that: (a) samples produced under tension by 
the Finkelmann procedure are intrinsically anisotropic and should show a small (temperature 
dependent) birefringence in the high temperature / phase, (b) for the I N transition via 
cooling there is a spinodal limit but for the N —* I transition via heating there is no soft mode 
at the standard spinodal temperature, (c) the N — > I transition is reminiscent of a martensitic 
transformation: nucleation of the I phase should occur in the form of platelets, making a well 
defined angle with the director. 



Introduction. - Nematic rubbers were first constructed in a single domain form by H. 
Finkelmann and coworkers 1 . They show a spectacular change in shape when they are 
switched from low temperatures (T < T c ) to high temperatures (T > T c ). Many properties 
resulting from the coupling between nematic order and elastic deformations |2] have been 
analyzed by M. Warner, E. M. Terentjev and coworkers We are concerned here mainly by 
the transition from the nematic phase N, to the high temperature phase /. 

Two striking facts should be mentioned here: (a) the transition is expected to be first order, 
but the plots of birefringence versus temperature T are continuous. 0] (b) the transition is 
very slow (minutes). [S] 

Our aim here is to discuss some effects of the nematic / elastic coupling on the phase 
transition. First, we discuss the anisotropy induced by the preparation method. We show 
that this preparation imposes a fixed external field coupled to the order parameter: under 
this field the plots of order versus temperature are expected to be continuous. Second, we 
analyze T-jumps: for a cooling jump we find a traditional spinodal transition. But for a 
heating jump we show that the elastic couplings tend to suppress the spinodal instability. 
Finally, we discuss the nucleation of an I phase inside a nematic single crystal, and show that 
it should take place in the form of platelets, with well defined geometrical conditions. 
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Anisotropy effects. - In the Finkelmann scheme, a very weak network is prepared first, 
and is put under a prescribed deformation e° zz along one axis. £Q Then a second reaction is 
started and the final nematic rubber is generated. The basic energy for this situation is of the 
form 

F = F (Q a p) - AQaflEfla + fj, (e aj g - £° a p) (spa - £° a ) (1) 

where summation over the repeated indices is implied. For simplicity, we restrict our attention 
to small deformations e a p. Fq is the standard Landau free energy of the nematic order, 
A describes the coupling between deformation and order, and the last term is the shear 
elastic energy (we consider only incompressible systems: Q a p and e a p are symmetric traceless 
tensors). The crucial point in Eq. Q is that the elastic energy (with a strong coefficient fi) 
is minimal in the original state (e = e ). We can shift the definition of deformations, writing 

e = e° + e (2) 

where e measures actual deformations from the high temperature / phase. Then the energy 
is (in contracted notation) 

F = F - AQsq - AQe + fie 2 (3) 

There is a constant field <jq = Aeo acting on the order parameter. This linear term in Q 
implies that, in the nominally isotropic I phase, we have a non-zero order Q. 
Writing explicitly for small Q 

F = A o {T)Q a0 Qp AQ (4) 

Via minimization in terms of e a p: 

e a p — ~2^Q a P (5) 
we first arrive at the classical renormalization of Aq(T) ~ ciq (T — To): 

F = a(T)Q aj3 Qp a - o-^Qpa (6) 

where 

a{T) = A Q (T) - = A (T) - a AT c (7) 
This minimized energy gives a finite high temperature order: 



cr. 



This describes a birefringence which is high near the transition point (where a is small) very 
much like the Kerr effect in standard nematics. This may be an explanation for the results of 
ref. g|. 

In the following, we employ a Landau expansion in terms of Q — > Q — Q 1 , and set the free 
energy to be zero at this "isotropic" phase, which now corresponds to Q = 0. 

Absence of spinodal instabilities. - Can we achieve the transition by soft phonon modes? 
Since our materials are essentially incompressible, we must investigate the transverse phonons. 
We assume rapid equilibration for the elastic degrees of freedom: we always use Eq. (JSJl. 

For a transverse phonon of wave vector q and displacement field u (with q ■ u = 0) 
the deformation is 

e a p = - (q a up + qpu a ) (9) 
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a) Consider first the instability of the disordered phase upon cooling (I — > TV transition) . 
Here the Landau free energy in Eq. © for small but nonzero amplitudes reduces to 

F = a(T)Q a0 Q Pa (10) 

and at the spinodal temperature T** the coefficient a(T**) vanishes. Note again that in this 
expansion, Q is actually the difference from Q 1 . 

Amplitudes of our transverse phonons are proportional to Q a p (Eq. (0) and the en- 
ergy required for the phonon generation vanishes at T — T**: the phonons are soft at this 
temperature. 

b) the situation is different for the N I transition upon heating. Here we start with 
a uniaxial nematic phase with nonzero components (Q X x 1 Qyy, Qzz) — i~S/2, —5/2, S); we 
can find a candidate for the spinodal temperature T* for a special value S = S*, satisfying 
the conditions 

OF d 2 F 

dS = dS^ = ( at S = S* and T = T*) (11) 
We employ the Landau expansion (which is already minimized for the elasticity) 

F — F — aQ a pQp a — bQ a pQp 1 Q ia + c (Q a pQpa) [Q^sQs^) (12) 

where b and c are positive (see Ch. 2 of ref. 0). 

We expand the free energy to second order around the spinodal point. At this point the 
coefficient d 2 F/dQ 2 zz — but the other curvature coefficients do not vanish. We find: 

% 2 r i 
F(Q) - F(Q*) = — [(Q xx - Q yy ) 2 + AQ%\ (13) 

The linear term in AQ a /3 = Q a p — Q* a p and the AQ^ Z term vanish as expected from Eq. 
i|llfl . In addition, there is no term proportional to Q 2 ZX and Q 2 y - they are ruled out by the 
rotational invariance around the director axis z. 

We can now investigate the transverse phonon modes (Fig. 

The first mode, with q and u in the x — z plane (Fig. ^l), induces a component propor- 
tional to Q xx (— q x Ux = esin0cos0 with e = qu), but not to Q vy . Thus, it contributes to 
the (Qxx ~ Qyy) 2 term in the free energy, which is positive (~ l^e 2 y ~ (esin6>) 2 ). A finite 
Qxx — Qyy implies a biaxial nematic, and this costs energy. 




Fig. 1 - (a) First transverse mode it = (it cos 6, 0, — usin 8). (b) Second transverse mode u = 
(0, u, 0). Note the difference from Fig. [5] the director vector is parallel to the z-axis. 



The second mode, with u normal to both q and director vector, induces an amplitude of 
Q xy (~ q x u y — esin#): again this induces biaxiality and requires a finite energy (~ ^e 2 y ~ 

(£sin6>) 2 ). 
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Fig. 2 Fig. 3 

Fig. 2 - (a) Habit plane separating nematic and isotropic regions. In two dimensions the "isotropic" 
network is represented by a square array, while the nematic network is represented by parallelograms. 
Along the habit plane the two networks match, (b) Detailed definitions of angles and axes. 

Fig. 3 - Platelet nucleus in the nematic phase accompanied by a large strain around the periphery. 



We conclude that both transverse modes are not soft for all nonzero angles (9 ^ 0) at 
T = T* . The = mode is soft at all temperatures, if F(Q) has full rotational invariance. 
But this mode cannot cause a change in the magnitude of the order parameter and does not 
catalyze the transition. 

Nucleation of a high temperature I phase. - We consider now a T-jump from a temper- 
ature just below the thermodynamic transition point T c towards a higher temperature T, in 
the region where spinodal instabilities are ruled out. 
A ) Choice of a habit plane 

We now investigate a possible plane boundary (habit plane) between / and N phase. 
A.l) The two-dimensional case (plane strain) 

This two-dimensional example (in the x — z plane) is shown on Fig. [21 where we impose 
compatibility between an / phase (represented symbolically by a square unit cell) and a 
uniaxial N phase (represented by parallelograms). The squares have been rotated by an 
angle <f> from the x-axis. The N phase has an elongation e zz along the z axis, which will be 
determined below. Note here that the director is not parallel to the z axis; the parallelograms 
are sheared with a strain 

e xz = ~(d x u z +d z u x ) = -tan.0 (14) 

Here and hereafter displacement fields are defined in the x — z frame shown in Fig. ^jp. There 
is also a deformation e xx = —e zz in the N phase (we assume incompressibility) and matching 
imposes 1 + e xx = cos <p. The two-dimensional deformation matrix is written as 



e xx e xz \ _ ( cos (j> - 1 \ tan < 



&xz &zz I V -^tan0 1 — cos( 



(15) 



Thus from Fig. |2h> the angle between the director axis and the habit plane is given by 

tan (9 + 0) = 1 ± e " ± Shl + = 2 = C ° S ± Sln (16) 
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Note that, for small (f>, the parallelogram is a lozenge, and the director is parallel to the long 
axis of the lozenge - at 7r/4 from the habit plane (Fig. [3] corresponds to this limit). 

The magnitude of cf> is fixed by the condition that the elastic energy of N state be equal to 
the equilibrium value f(Ti) in the initial nematic phase. Under the incompressibility condition, 
we have /(!») = fi(e^. x + e 2 zz + 2e 2 z ) = 2(i(e 2 . x + e 2 xz ). Thus, 4> is determined by 

/(H)/(2/i) = (1 - cos 0) 2 + 1 (tan <f>) 2 (17) 

For this two-dimensional case, we achieve a perfectly isotropic state / on one side and the 
standard N state of the other side. However, as we shall see, this perfect matching is not 
possible in 3D (exactly as in similar martensitic transformations [7|). 
A. 2) The three-dimensional case 

We can try a similar construction for the three dimensional case (see Fig. [21 again). We 
take the habit plane perpendicular to the x — z plane (the director axis stays in the x — z 
plane). In the rectangular coordinate composed of the director axis (n), y axis, and a axis 
(perpendicular to the others, see Fig. [2J), the non-zero strain components are e n (T i ) = ^ 
(in the director axis), and e aa — e yy — — e,/2. The elastic energy associated with this is given 
by fi = f A* e l- I n the (x, y, z) coordinate the strain tensor is given as 

e xx e xy e xz \ / cos0-l ^tan< 



e V x e-yy e yz = -e%/2 I (18) 

e zx e zy e zz ) \ |tan^ 1 - cos 4> + e,/2 

and the elastic energy is given by = g xx + e^ y + e 2 zz -\-2e 2 zx . Thus, we can fix <f> by requiring 
that the energy is the same in the (n, y, a) and (x, y, z) coordinates: 

^e 2 = (cos $ - i f + (1 - cos <t> + e,/2) 2 + - (tan <f>f (19) 

In addition, we have 

2 — cos <b + e» 12 + sin 6 
tan (6 + 6 ) = Y 20) 

COS(p 

We note that Eqs. I|19|) and (|20|) predict that the habit angle approaches 7r/4 again in the 
small deformation limit (e^, <f> — * 0). 

In this three-dimensional case, the N phase has a third component of strain e yy = — e^/2 
and this can not be eliminated in the / phase. From the incompressibility condition in the 
I phase, we must also have a strain epp = —e yy ; inside the nucleus, there must be non-zero 
strain components {e yy ,ef3p) = (— e.;/2, ej/2), which results in a residual elastic energy (per 
unit volume) 

Ei = Me 2 /2. (21) 

Thus we reach a state which is not totally isotropic 
B) Homogeneous nucleation 

Again we think of a T-jump from an initial temperature Tj just below T c (we take Ti = T c 
in practice for simplicity) up to a final temperature T lying about 10 degrees above T c . Our 
assumed shape for the nucleus is a platelet corresponding to the habit plane orientation, with 
a thickness h and dimension L (Fig. |3J)- From now on we construct only crude estimates, 
ignoring all numerical coefficients: for instance the platelet volume is taken to be L h. 

The platelet energy is of the form 

/ = -L 2 h ■ A + E!L 2 h + 1 L 2 + E p (22) 
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Here A is the gain in free energy (per unit volume) obtained by switching from N to I at the 
temperature T. Thermodynamics imposes 

A = (T - T c ) AS (23) 

where AS ~ aoQ 2 is the entropy jump at the transition temperature T c . 

The second term Ej in Eq. I)22|) is due to the residual stress in the / phase [The principle 
axes of the deformation here are the y and (5 axes and another axes normal to them, i.e. the 
non-zero elements of the diagonalized deformation matrix are {e vv ,epp) — (— ej/2, e%/2)] and 
is only a fraction of the initial elastic energy. The original elastic energy is §/xe 2 while we 
have here £t(e^ + e|g) = |/xe?. 

Thus we lump the elastic correction into A, writing 

A = A - i/xe 2 = (T — T m ) AS (24) 

where we have introduced T m in analogy with Eq. ilL'-!l) . T m is the minimal final temperature 
for nucleation. From a scaling point of view, we should have 

A T-T m 

T4 " SST (25> 

where AT C corresponds to an increase in the N — I transition temperature due to the coupling 
(See Eq. 0). This can be explicitly shown by using the relation /ze 2 ~ A 2 Q 2 /fi and AS* ~ 
aoQ 2 . We assume from now on that A is positive but small. 

The third term in Eq. I)22|) is a surface term, with a certain interfacial energy 7, which 
might be estimated by a Ginzburg - Landau type free energy (including the spatial dependence 
of the order parameter) plus elastic tensors. Because of these latter terms 7 is positive even 
when T > T* as indicated by the previous arguments on the absence of the spinodal instability. 

The last term in Eq. (1221) is another elastic term, due to distortions at the periphery 
of our platelet. Here, there is no matching at all (the matching was achieved only at the 
habit plane boundary) , we have large deformations taking place in a toroidal region near the 
periphery (Fig. The length scale relevant for this extra peripheral deformation and the 
size of deformation are both around h and, because the deformation fields always satisfies the 
Laplace equation, this deformation dies out only at a distance ~ h from the torus: the strain 
e p (at least of the order of e^) is stored around the periphery of the volume ~ Lh 2 (Fig. (SJ: 

E p ~ Lh 2 A p (26) 

with 

A p = [ie 2 p (27) 
We now write Eq. 1221) in the following form: 

F = -AL 2 h + 7L 2 + Lh 2 A p (28) 

By minimizing the above energy with respect to h, we have 

ho - (29) 

F = -^AL 3 + 1 L 2 (30) 
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Thus the assumption L 3> h corresponds to A < A p : if the peripheral energy A p is large 
compared with the bulk gain A, the platelet tends to become thin. 

The energy expression in Eq. (|3U[) is in contrast with the classical case of a spherical 
nucleus without any elastic contribution: F = —A • L 3 + jL 2 . For A p > A, the bulk term for 
the disk droplet is quite small compared with the surface term. If we look for the maximum 
of Eq. (|3L)[1 with respect to L, we find an energy barrier U = Fq(L*) for the formation of a 
critical droplet: 

v _ 7 8A P 
A 3A 

3 V 3A 2 J 

where U c ~ 7 3 /A 2 is the standard barrier. In our regime U 3> U c the nucleation rate can be 
thus dramatically suppressed. 

Conclusions. - (1) Eq. (|32[) shows clearly that homogeneous nucleation is prohibited in 
nematic rubbers - just as it is in martensites. Clearly, we need heterogeneities to nucleate, 
and we do expect them in networks. (2) Even when dealing with heterogeneous nucleation, it 
is reasonable to think that platelets (or needles) close to the habit angle will be preferentially 
generated. This would possibly be observed by electron microscopy or by force microscopy on 
the outer surface. (3) From a practical point of view, it should be beneficial to accelerate the 
commutation process - we can think of at least two ways: (a) adding colloidal platelets (clay 
?) at the correct angle (possibly orienting them by fields during the synthesis), (b) shearing 
the sample with a shear plane near the habit angle. 
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